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We have constructed the lowest few eigenvectors of the staggered Dirac operator on 5(7(3) gauge configurations, 
both quenched and dynamical. We use these modes to study the topological charge and to construct approximate 
hadronic correlators. 



1. Introduction 

In a recent study jlj] we successfully calcu- 
lated the mass of the rj' meson using staggered 
fermions. According to conventional lore, the r]' 
receives a large portion of its mass from instanton 
effects. At the fermionic level, this would mean 
receiving contributions from zero modes of the 
Dirac operator. These zero modes come with a 
definite chirality in the continuum and are asso- 
ciated with the topological charge (Q) of a gauge 
configuration via the index theorem: 

Q = n+ - n_ (1) 

where n + and n_ are respectively the number of 
right and left-handed zero modes. At finite lattice 
spacing, 75 and p do not exactly anti-commute, 
so there are no exact chiral zero modes and nor 
is the topological charge of gauge configurations 
well defined. Nevertheless, one might expect to 
find a few low lying eigenmodes carrying chiral 
charge. Smit and Vink j^] verified this sort of 
behavior for 17(1) gauge theory in 2 dimensions. 

In this work, we set out to study the extent 
to which the lattice "zero modes" reproduce the 
continuum picture for QCD. Accordingly, we have 
constructed the lowest few eigenmodes of stag- 
gered p on both quenched and dynamical SI/ (3) 
gauge configurations (83 samples each) of size 
16 3 x 32 corresponding to a lattice spacing of 
O.lfm. We have used the subspace iterations 
method Q to compute the eigenmodes. The 
eigenvalues that we obtain have relative errors in 
the range 10~ 7 for the smallest through 1CP 4 for 
the highest. 



2. Topological charge 

The lattice analog of the U(l) axial anomaly 
equation gives 

Q = rntr(T 5 G( X ,y))^rnj: ^^ (2) 

where G{x,y) is the fermion propagator. The 
mode by mode contribution to u\T5U\ on all 
the dynamical configurations is shown in Fig- 
ure |l|. It appears to diminish as A increases which 
can be verified quantitatively by calculating the 
width, shown in the bottom part of figure [l]. 
This behavior is consistent with the continuum 
where u^r 5 it>, = for all the non zero eigenval- 
ues. Figure || shows a similar set of plots for the 
quenched ensemble. Although both (T^) and its 
width show qualitatively similar behavior as be- 
fore, they seem to drop more slowly with increas- 
ing A. The solid line present in all the figures 
shown so far represents the minimum value of the 
32nd eigenvalue in the ensemble, ie. the line be- 
low which we are sure to have found 100% of the 
modes. The sum in eqn(||) for r = 2, 4, 8, 24 
and 32 modes at ma = 0.01 on a typical dy- 
namical configuration is shown in figure [|. The 
plateau in Q as the number of modes increases 
implies that the contribution of the higher modes 
to the sum in eqn|| is becoming negligibly small. 
Although, the plot shown here is for one con- 
figuration, a similar behavior is seen on all the 
configurations, quenched and dynamical. In Fig- 
ure ||, the isolated point represented by a square 
is the value of Q obtained from pseudofermions. 
The "true" value obtained by averaging over 82 
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Figure 1. Upper:(r5) on the dynamical ensemble. 
Lower: Fluctuations in (T5). 



copies of noise lies within the band marked by 
the dashed lines. Since the plateau in Q lies 
within this band it is a very good indication that 
these lowest modes are the topological zero modes 
which make the dominant contribution to the sum 
in eqn ||. We have compared Q obtained from 
pseudofermions and modes on all the dynamical 
configurations and find that both agree within er- 
rors. On the quenched ensemble, the difference 
between the two estimators appears to be statisti- 
cally significant. This simply reinforces the infer- 
ence reached by comparing Figures |l| and ^| that 
more than 32 lowest eigenvalues are required to 
identify an upper limit for the number of shifted 
zero modes on the quenched configurations. 

3. Hadron Correlators 

Calculation of hadron correlators on the lat- 
tice normally involves inverting p + m to obtain 
the quark propagator. When resolved in terms of 
the eigenvalues and eigenvectors of p, on a fixed 
background gauge field, it takes the form 
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Figure 2. Above: (T5) on the quenched ensemble. 
Below fluctuations in (r 5 ). 
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Figure 3. Topological charge on a typical dynam- 
ical configuration 

Comparing the full propagator obtained using 
standard inversion procedures like CG, with that 
constructed out of the sum above with the avail- 
able eigenmodes gives an insight into the role of 
the eigenmodes in determining hadronic quanti- 
ties. 

The rf meson, being a flavor singlet, has 
both connected and disconnected contributions 
to its correlator. The discussion in section 2 
suggests an extension of the analysis to the 
Tj disconnected correlator which is given by 
(Tr(r 5 G(x, x)) Tr(r 5 G(y, y))). We calculate this 
quantity using pseudofermions and with the avail- 
able 32 modes. The ratio R(t) of disconnected to 
the connected correlator for ma = 0.01 is shown 
in fig H For the connected correlator, we use the 
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full result available from our previous study |Q. 
On the dynamical configurations, R(t) is fit to the 
form § R(t) =fy[l- exp(-(ffty - m s )t)]. For 
staggered fermions, N v , the number of valence 
fcrmions equals 4 while m\ = (4m^- + 3m^. + 
mZ)/8. It is evident from the figure that R(t) 
data with 32 modes agree within error with the 
full result. We find similarly good agreement 
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Figure 4. Pion-?/ splitting. 

with only 24 modes though 16 appears to be not 
quite enough. The pattern is the same at the 
heavier quark masses, 0.02 and 0.03. 

We note that this behavior is not automatic 
as 32 modes are insufficient to yield a good ap- 
proximation to the true quark propagator. We 
find, for example, that the pion is off by a large 
factor. Based on our knowledge of the staggered 
eigenvalue spectrum and on recent work done by 
Negele et al. Q with Wilson fermions, one might 
expect that 100-200 modes would be needed to 
reproduce the pion propagator. Thus the r]' is 
special in exhibiting an extraordinary sensitivity 
to the lowest few modes. 

The quark propagator, constructed out of the 
available 32 modes, can be used as a starting seed 
for calculating (Jj) + m) _1 using CG. The benefit 
of using such approximate solutions as a starting 
iterate in the CG procedure is illustrated in Fig- 
ure ^. The number of CG iterations required for 
the residual to converge to a given tolerance drops 
by more than 50% at ma — 0.01. At higher quark 
masses, this effect is not so pronounced owing to 
the fact that the lowest 32 eigenvalues are smaller 
than the quark mass itself. These conclusions 
hold only for delta function sources or gauge fixed 
wall sources built out of delta functions. Quark 
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Figure 5. Comparison of CG convergence. 

propagators from noisy sources are dominated by 
short distance modes and hence there is virtually 
no improvement in the CG convergence rate when 
the starting seed is built out of long wavelength 
modes. 

4. Conclusions 

We find that it is possible to identify the 
shifted zero modes on the lattice since they 
make the dominant contribution to uxT^ux and 
tr(TsG(x,y)). On our dynamical configurations, 
it seems that the number of topological zero 
modes must be less than 32 while on the quenched 
configurations the data indicates that there are 
more than this number. For most hadronic quan- 
tities, the approximate propagator constructed 
out of these topological zero modes is not accu- 
rate enough to be useful. However, a remarkable 
result of our study is that we find that the rf fla- 
vor singlet interaction on the dynamical config- 
urations is dominated entirely by the quasi-zero 
modes of p. 
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